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Systems under holonomic constraints are classified within the generalized Hamiltonian framework 
as second-class constraints systems. We show that each system of point particles with holonomic 
constraints has a hidden gauge symmetry which allows to quantize it in the original phase space 
as a first-class constraints system. The proposed method is illustrated with quantization of a point 
particle moving on an n — 1-dimensional sphere S" 1-1 as well as its field theory analog the 0(n) 
nonlinear sigma model. 
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A systematic study of the generalized Hamiltonian dy- 
namics was made by Dirac |l| who classified the con- 
straint structure and developed the corresponding opera- 
tor quantization schemes of the constrained Hamiltonian 
systems. 

Gauge theories appear under first-class constraints. 
The Dirac's theory of first-class constraints was com- 
bined with the Feynman path integral method HQ]. A 
general approach to quantization of first-class constraints 
systems was developed byFradkin and his collaborators 
0, 0, it is reviewed in [||. An alternative symplectic 
scheme has been proposed by Faddeev and Jackiw 0, @ • 

Among physical systems under second-class con- 
straints are anomalous gauge theories 

H IIS E3, the 

nonlinear chiral model [l2j, many-body systems involv- 
ing collective and independent-particle degrees of free- 
dom [l3L IT3. ITU - The Dirac's quantization method of 
such systems consists in constructing operators replicat- 
ing the Dirac bracket for canonical variables and taking 
constraints to be operator equations. Batalin and Frad- 
kin ^(| proposed a quantization scheme which converts 
constraints to first class by introducing new canonical 
variables. The problem reduces thereby to quantization 
of a first-class constraints system in an enlarged phase 
space. This method has found many useful applications. 

The Hamilton- Jacobi procedure is also used for quan- 
tization of constrained systems |f?]> Il8l Il9| . 

We show that each system of point particles with holo- 
nomic constraints has a hidden gauge symmetry which 
allows to perform its quantization without introducing 
new canonical variables. A mathematical pendulum on 
an n — 1 -dimensional sphere as an example of a 

mechanical system and the 0(n) nonlinear sigma model 
as an example of a field theory are quantized using the 
proposed method. 

Let C be Lag rangian of point particles with holonomic 
constraints 



where Xa are Lagrange multipliers, xa = XA{4> a ) = 0, 
with A — f ..771, are constraints in the coordinate space 
(<j> a ) with a = I..71, and U — U((f> a ) is a potential energy 
term. 

A more general class of systems with kinetic energy 
T = ^^2a i j{q)q l q : > can be converted to |QJ by an iso- 
metric embedding into a euclidean space of a higher di- 
mension plj . 

The vectors e a A = dxA/d(j) a constitute a local basis in 
a manifold orthogonal to the constraint surface xa = 0. 
By linear transformations of xa with 0-dependent coeffi- 
cients, one can ensure the global orthogonality conditions 

e a A e% -Sab M 

A straightforward follow-up to the Dirac's scheme 
leads to the following results: Let ir a and tta be canon- 
ical conjugated momenta corresponding to the canon- 
ical coordinates (f> a and the Lagrange multipliers A^, 
respectively. The primary Hamiltonian is given by 



H p = 



l 7r" + U — \aXa- The primary constraint 
functions are Qa — ^ A- The canonical Hamiltonian 
is H c = H p + uatta where ua are undetermined func- 
tions of time. The holonomic constraints \a — appear 
again within the generalized Hamiltonian framework as 
secondary constraints: xa = {Qa,H c } ~ 0. The con- 
straints of the next generation look like VLa = {xa, H c } = 
TT a e A k 0, while the third generation is va = {^A, H c } = 
T: a T^d 2 XA/d(j) a d(l) 13 - e'\dU/d(t) a + e a A e%X B « 0. Equa- 
tions va = allow to fix A^4 in terms of <jf and -K a , no 
new constraints then appear. 

The phase space can be reduced by eliminating the m 
canonical pairs (Xa,t^a) from the Hamiltonian and the 
Poisson bracket. We solve equations Qa = with respect 
to tta and va — with respect to Xa and substitute 
solutions into H c . The result is 



d 2 XA 



dU 
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U + XaXa 



(1) 



K = ^ a + UH^^^-e a A^)(e a A e B )- 1 XB. 

(2) 

The remaining constraint functions xa and Qa sat- 
isfy {xa,Xb} = 0, {xa,^b} = e%e%, and {MaMb} = 



2 



deP A /d<j> a e B 



P. 



7r" — Oe' B / Ocfre A -K a . These are second-class 
constraints. The H' c is first class. 

Two sets of the constraint functions are equivalent if 
they describe the same constraint surface. Two hamil- 
tonian functions are equivalent i.e. provide the same 
hamiltonian dynamics if they and their first derivatives 
coincide on the constraint surface. E.g. the replace- 
ment (e^e^) -1 — > Sab in Eq.JSJ results to an equiv- 
alent Hamiltonian, since the variance is of the second 
order in the constraint functions. The replacements 

also, etc. 

The replacement Q, A — > Q' A = Cl A — de A /d(j> a e B ir a XB 
gives an equivalent set of the constraint functions Q a = 
(xa, Q'a) with the Poisson bracket 



— n a ~ e A &,A provide a possible modification 



{Qa, Qb} = lab + C- a b Qc 



(3) 



where Tab = 1 m +A m+B = 0, 1a m +B = —l m +B a — 
Sab, so that J ab = —l ba and T ab Z hc = —S ac . In what fol- 
lows, tensors with the upper and lower indices are related 
by % = l ab T b . 

The constraint and hamiltonian functions can be ex- 
trapolated from the constraint surface Q a = into the 
whole phase space in various ways. We wish to establish 
the normal form of the constraint and hamiltonian func- 
tions, Q a and Ti., for which the relations {Q a ,Qb} = ^ab 
and {Q a , Ti} = are valid in the strong sense, at least 
locally near the constraint surface. 

Let us set Q [ a 1] = Q a , Ti [1] = Ti' The relations 



{Qa k] ,Q l b k] } = lab^ 

{Q[ k \H^} = Hp- 



n ci 
^ab 



bk G, 



■Ckn[k] n[k] 
=>ci ■■■Vc k i 

Q [k] 



(4) 
(5) 



are valid for k = 1. If Eqs.0} and @ arc fulfilled for k > 
1, the equivalent constraints and hamiltonian functions 
of order k + 1 can be constructed: 

The Jacobi identities and symmetry of the structure 
functions C^ 1 ' " Ck and TZ bl "° k with respect to permuta- 
tions of the last k indices give 



Haby...b k — TZb ia ...b k ~ 0. (7) 



again. At each step, the Poisson bracket relations get 
closer to the normal form. In the limit k — ► +oo, we ob- 
tain {Q a , Qb} — lab and {Q a , Ti} — 0. A similar theorem 
on the normal form exists for systems under first-class 
constraints 6] . 

The xa depend on the <j> a only, Q! A ~ 7r Q , and Ti' c = 
T + V where T ~ n a w^ and V ~ 1. These properties 
propagate throughout the above algorithm, being thus 
inherent for Q a and Ti: 

Indeed, let I,J,K.. = l..m and Q, R,S.. = m + 1..2m. 
Let at an k-th step Qi be a function of ip a and Qq be 
linear in ir a . By inspection of Eq.(@J one can conclude 
that Cj Qci - Ck = for all c t , cf Cl " Ch ^ only if c % = 
K, L, M... The coefficient Cg Cl "' c * can be expanded into 
zero- and first-order terms in 7r Q . The zeroth term does 
not vanish if a — J,K,L.. at i = l..k except for i = I 
with ci — R, S, T... The first-order term does not vanish 
if Ci = J,K,L.. at i = l..k. Eq.JBJ implies that Qi ~ 1 
and Qq ~ 7r a at step fe + 1. 

Eq.JHJ) is valid at k = 1 separately for T and V. Let 
~ 7r Q 7r' 3 . By inspection of Eq.QJ one finds that 
structure functions of the kinetic energy T obey TZj 



J...K 



QJ...K 



QRJ...K 



o, n. 



J...K 



-r> _RJ...if 



1, and n Q HSTJ - K = 0. 



1, ft* 

7r a , n Q HSJ - K 

Eq.© implies that the next-order kinetic part of the 
Hamiltonian remains quadratic with respect to ir a . Let 
further ~ 1. Eq.QJ tells that structure functions 
of the potential energy V satisfy TZq J,k ~ 1, the other 
functions vanish. In virtue of Eq. © , the next-order V is 
again a potential. 

The above algorithm assures Q a ~ Q a , 'H'c ~ H, 
{/,^a} « {f,Qa}, and {/,«;} «_{/,«} for arbitrary 
functions /. In the basis Q a = (xa, &a)i E a — dxA/d(f) a . 
Since Qa is linear in the momenta, it can be written as 
Ha = Tr a F A . The normal form of the constraints implies 
E%Fg = Sab and 



F^dFg 



(10) 



E a A 



F%. 



In addition, we have e A 

The constraint functions Q a generate transformations 
of the canonical variables 



The next-order constraint and hamiltonian functions are 
given by 

g[k+l] = e W + _l_ Co Coc 1 ...c^M W^W ) (8) 



[fe+1] 



n 



ik] 



L -K b ^- bk Q k }Q k] ...Q k }. (9) 
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They define the same constraint surface and the same 
hamiltonian flow on the constraint surface as one step 
before. Using Eqs.Q) and (JSJ, one can calculate the 
next-order structure functions C f 1 '" Cfc+1 and TZ a x " h+1 
and repeat the procedure. If structure functions vanish, 
we shift the k by one unit and check Eqs.(0} and |J5J 



sr = {<f> a ,XA}e A = o, (ii) 

Sir a = {n a ,x A }6 A = -E%9 A , (12) 
= {<t> a ,nA}r,A = F2 VA , (13) 
At: = {7r a ,Cl A }VA^-^dF^/drVA, (14) 

with respect to which Ti. is invariant, with 8a and t\a 
being infinitesimal parameters. 

Lagrangian C(4> a , ir a , 4> a ) = n a (j) a — H calculated for 
<p a = {cj)"^} is invariant under l(TT|) - (fT^f> up to a 
total differential of some function. Suppose we have 
solved the </> Q = {>", H} for the n a and got £(cj) a ,<j) a ) = 
jC(4> a , ir a {(j)^, cf)P), 4> a ). The transformation properties 
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the (jf are as follows 

= {{</> a ,H},XA}6 A = 0, 
A0« = {{(/> a ,H},n A }i lA = dFZ 



A/3 



(15) 
VA- (16) 



Lagrangian C(4> a ,(j> a ) as a function in the tangent bun- 
dle {4> a , 4> a ) is invariant with respect to the m-parameter 
family of the global transformations 113|) and (|16fl . Equa- 
tions {CIat'H} = are equivalent to the following equa- 
tions 



d<j) a A 9 







(17) 



which indicate explicitly 23j the invariance of £ against 
the transformations (|13[) and (|16|) . 

The constraint functions xa do not induce transfor- 
mations in the tangent bundle (cf> a ,(j) a ). This is a spe- 
cific feature of holonomic systems whose constraints do 
not depend on velocities. Equations {xa,W} = im- 
ply, however, E^fi* = 0. Consequently, C((f> a , <j) a ) is 
defined on the constraint surface E\<j> a = 0. It can 
be extrapolated from the constraint surface according to 



C'(4> a ,(p a ) = C((f> a ,A a P(pP) where A afj = 8 af) - F%E% 
A Q/3 acts as a projection operator. £' is invariant under 
an m-parameter family of transformations 



8'<j> a = 0, 
6'<p a = F%§ A 



(18) 
(19) 



where d A are infinitesimal parameters. Equations Q A 



n apa _ q are e q U i va i en t to the following equations 



dC 



F% = 



(20) 



which show the invariance of C under transformations 
(|18fl and l|19|) . Lagrangian C, as a Legendre transform 
of H, depends thereby on the tangent velocities A Q ' 3 </>' 3 
automatically, i.e., C'{<j) a , <j) a ) = £(</>",</>"). 

The invariance of Lagrangian C under the 2m- 
parameter family of the global transformations (| 1 31) and 
(|16[) . I|18l) and l|19l) provides the invariance under an m- 
parameter family of local transformations of trajectories 
4> a {t). Let T) A be time dependent functions and i) A = f} A . 
The transformation law in the tangent bundle (0 Q ,0 Q ) 
reads 



5'<j) a + Acj) a 



(21) 

VA = | (F2va) (22) 



The last equation suggests that this transformation acts 
on the trajectories 4> a (t). Lagrangian C is therefore in- 
variant with respect to the local transformations. The 
variations of 4> a (t) along the vectors F% correspond to 
the m gauge degrees of freedom. The constraints xa — 



play a role of the gauge-fixing conditions. The gauge- 
invariant observables are projections of trajectories </> a (i) 
onto the constraint surface xa = 0. Ea. 1)10(1 ensures that 
the symmetry group is abelian. 

If C would be chosen as the starting Lagrangian to 
pass over to the Hamiltonian framework, one could get 
TL as primary Hamiltonian and Qa = Tt a F2 — as pri- 
mary first-class constraints by virtue of Eq. I|20|l . The con- 
straints xa = can be selected as a set of the admissible 
gauge-fixing conditions. Starting from Lagrangian 
we arrive at the set of second-class constraints xa = 
and fl A = 0. Starting from C, we obtain the same set as 
gauge-fixing conditions xa = and first-class constraints 
&a = associated with the gauge symmetry. It shows 
that classification of constraints in the holonomic systems 
is a matter of convention. 

As Lagrangian 0J) is equivalent to its gauged counter- 
part C(4> a , (j) a ) under the gauge-fixing conditions xa = 0. 
one can apply one of the quantization schemes specific 
for gauge theories. 

The system is quantized by the algebra mapping 
(0 Q ,7r a ) ->■ (0 Q ,7r a ) and {,} -> Consequently, to 

any symmetrized function in the phase space variables 
one may associate an operator function. The function 
is symmetrized in such a way that quantal image is a 
hermitian operator. 

The wave functions obey the Dirac's supplementary 
condition 



n A ^ = o. 



(23) 



In the quasiclassical limit, gradients of the wave functions 
along the vectors F% vanish. 

The path integral for the evolution operator becomes 



Z = 



(24) 



J Y[d n 0d n nl[S(x A )5(fi A )exp{i J dt^ a ^ a ~fL) 

Eqs.J2BJ and 124|) solve the quantization problem for sys- 
tems of point particles with holonomic constraints. 

Given the initial state wave function \t(0) satisfying 
(|23p. the final wave function ^>(t) can be found applying 

the evolution operator l|24|l on 'J'(O). Since fl A commute 
with the Hamiltonian, the final state wave function obey 
Eqs.HSJ. 

It is desirable that the quantization procedure does not 
destroy the classical symmetries which results in having 
the supplementary condition satisfied by the state 
*&(t) for any value of t. This feature can apparently be 
violated either due to complex terms entering TL or by 
approximations adopted for treating the operator eigen- 
values. The supplementary condition (|23[) assures, how- 
ever, that the symmetry is preserved in average. As we 
shall see latter on, the quantization procedure preserves 
the gauge symmetries explicitly in the applications con- 
sidered. 
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The supplementary condition i|23|) does not comprise 
the constraints xa — 0- In gauge theories, the evolution 
operator is independent of the gauge-fixing conditions 
0| . We can insert detjx^^B} = 1 into the integrand 
of Ea. l|24|l to bring the measure into the form identical 
to the one prescribed by gauge theories. In the path in- 
tegral, the xa can then be replaced with arbitrary func- 
tions. 

It is remarkable that physical observables are deter- 
mined only by half of the number of second-class con- 
straints. 

Let us illustrate the method proposed with quantiza- 
tion of a mathematical pendulum on an n— 1-dimensional 
sphere 5™ -1 . 

Lagrangian is given by Eq. Q with m = 1 and x = In 4>- 
Using the Dirac scheme, we obtain the primary Hamil- 
tonian 7i p = ^ir 2 — X\n(f> and the primary constraint 
Ga = 7T\ w 0. The canonical Hamiltonian becomes 
7~Lc = "Hp + U7r \ where u is an undetermined function 
of time. The secondary constraints Ga+i = {G a ,T~tc} 
can be found: Gi = Incf), Qi = (jyir/cf) 2 , G 3 = n 2 /(j) 2 — 
2(0tt) 2 /(/) 4 + \/<p 2 . The last constraint G 3 = allows 
to fix A in terms of 4> a and 7r Q , no new constraints 
then appear. The canonical pair (A, ir\) is eliminated 
by solving equations Go = and G3 = 0. The result is 
H' c = ±Tr 2 + (n 2 -2((f>Tr) 2 /<p 2 )\n(f>. There remain two con- 
straint functions Gi and G2, such that {Gi,G2} = 1/0 2 - 
These are second-class constraints. The Hamiltonian TC' C 
is first class. 

In order to split the constraints into a gauge-fixing con- 
dition and a gauge generator, we should construct Ga and 
Ti. starting from Ga and Ti.' c . The same constraint surface 
is described by the functions x — In <f> and fi = </>7r satis- 
fying {x, f2} = 1, so one can set Gi = ln</> and G2 — 4 111 '■ 

The Hamiltonian TL can be found to be 

TL = i</» 2 A Q/3 7r Q 7 r' 3 (25) 

where A af3 = 8 af> - (j) a ^ / <t> 2 . For n = 3, H = \L a L a 

where L a is the orbital momentum. The difference TL — 
TL' C is of the second order in the constraint functions. 
This is sufficient to have identical hamiltonian flows on 
the constraint surface. 

fl is stable with respect to time evolution {fl,TL} = 
0. The relations Acj) a = {<j> a ,n} = <j) a and Att q = 
{7r a ,17} = — ix a show that ft generates dilatations of 4> a 
and contractions of ir a . The x is identically in involution 
with the Hamiltonian {%, TL} = 0. The Poisson bracket 
relations {4> a , x} = and {tt , x} — —cj) a define shifts of 
the longitudinal component of the canonical momenta. 

The gauged Lagrangian is given by 

i = \^{m a ^i<t> 2 . (26) 

It is invariant under the local dilatations of <f> a . In the 
context of the dynamics defined by C with no constraints 



imposed, <j> is an arbitrary function of time. It can al- 
ways be selected to fulfill the constraint x = or some 
other admissible constraint. The Euler-Lagrange equa- 
tions have the form 

A°*»|^/0) = O. (27) 

The point particle moves thereby without tangent accel- 
eration. Eas. l|27|l are not affected by constraints. 

C depends on the spherical coordinates 4> a /4> which lie 
on an n — 1-dimensional sphere of a unit radius. Eas. l|27|) 
for unconditional extremals in the coordinate space <f> a 
coincide therefore with the D'Alembert-Lagrange equa- 
tions for extremals in the spherical coordinates space, 
considered as independent Lagrange variables, under the 
condition that they lie on an n — 1-dimensional sphere of 
a unit radius. 

The constraint x — can therefore be treated as a 
gauge-fixing condition, the function <j> as a gauge degree 
of freedom, the ratios <j) a /<p as gauge invariant observ- 
ables. The equations of motion l|27|) are formulated in 
terms of the gauge invariant observables. The extremals 
of the action functionals J Cdt and J Cdt under the con- 
straint x — coincide. 

The quantum hermitian Hamiltonian has the form 

TL = i^A^Tr'VA^TrT. The vector i0A a %^ gives the 
angular part of the gradient operator. Although it is 

not conspicuous, TL does not depend on the radial co- 
ordinate cj>. The constraint operator can be defined as 

fl = (<f) a Tr a + 7r Q Q )/2. It acts only on the radial com- 
ponent of cf) a , so the relation = holds. The sup- 
plementary condition (J2HJ implies ^ = 4>~ n / 2 ^i(<fi a /(f>). 
The important information about the system is contained 
in $ \{<j> a I '(f>) only. 

The integral over the canonical momenta Eq. l|24|l can 
be evaluated explicitly to give 

Z= [ I] V(9x/9r) 2 5(x)^exp{i / dtc). (28) 



Lagrange's measure y (dx/ dcf) 01 ) 2 5(x)d n <fr coincides with 
the volume element of the 5™~ 1 sphere. It can be rewrit- 
ten as an invariant volume of the configuration space e.g. 
in terms of the angular variables (<pi...<^ n _i) with the 
help of the induced metric tensor. Lagrange's measure is 
invariant under the group 0(n). It remains the same for 
all functions x vanishing for 0=1. 

The 0(n) nonlinear sigma model represents the field 
theory analog of the spherical n — 1-dimensional pendu- 
lum. The n — 4 case corresponds to the chiral nonlin- 
ear sigma model which is invariant under the SU (2)^ ® 
SU(2)r ~ 50(4) group. Its quantization does not in- 
volve essentially new ingredients. The coordinates <j> a 
become functions of x p — (t, x, y, z), the time derivatives 
<j) a in Eqs.QJ and (|26|l should be replaced with <9 M a . 
The wave functionals depend on the spherical coordinates 
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4> a (x) I <fi(x) , the path integral has the form of Eq.J33J) 
with dt replaced by d 4 x. 

The quantized gauged pendulum and 0(n) nonlinear 
sigma model preserve the hidden gauge symmetries. Con- 
sequently, the wave functions \& satisfy the Dirac's sup- 
plementary condition at any time, i.e., the systems evolve 
always on the constraint surfaces. 

We showed that systems of point particles with holo- 
nomic constraints have hidden abelian gauge symmetries. 
Their physical origin is the following: We allow virtual 
displacements of the particles from the constraint sur- 
face and treat such displacements as unphysical gauge 
degrees of freedom. The holonomic constraints turn 
thereby into gauge-fixing conditions. The physical vari- 
ables are fully specified by projecting the particle coor- 
dinates onto the constraint surface. In the considered 
examples, the projected coordinates are the spherical co- 
ordinates <j) a '/</>. The existence of the equivalent gauge 
models allows us to quantize the original holonomic sys- 



tems in their own phase space using standard methods 
of quantization known from the gauge theories. 

After submitting this work, our attention has been 
called to papers Hi |H El El IH . In Ref. 0, 
a mathematical pendulum is quantized using the Dirac 
approach. The 0(n) non- linear sigma model is discussed 
along the line in Ref. pfij . In Ref. pfij , hidden gauge 
symmetries of second class constraints systems are dis- 
cussed and a quantization method in the original phase 
space and an extended configuration space is proposed. 
The holonomic systems, as we showed, admit the quan- 
tization in the original phase and configuration spaces. 
The methods 0, |2g are applyed for quantization of var- 
ious physical systems and compared in Refs. [27l EsL l2flj . 
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